Parent phases of doped Mott insulators: a Cluster DMFT study 
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We investigate the insulating phases of a frustrated Hubbard model in its strong coupling regime at 
half-filling. We pay special attention to all the symmetry breaking instabilities that can be described 
by Dynamical Mean Field Theory (DMFT) on a square plaquette. We identify several mean-field 
solutions, two Neel states breaking the SU{2) symmetry, a dimer phase and a Mott insulator which 
doesn't break any symmetry. The singlet to singlet fiuctuations soften dramatically in the latter 
phase, giving rise to dimerization fiuctuations as well as chiral fiuctuations that are both low-lying. 
We present a simple picture of the different DMFT states and their evolution with frustration. 



Frustrated two dimensional Mott insulators are inter- 
esting in their own right, and are realized in Helium films 
as well as in transition metal and organic compounds. 
They received enormous attention, following the discov- 
ery of superconductivity in copper oxide based materials 
and the suggestion [ ] that the phenomenon of high tem- 
perature superconductivity is connected to the doping of 
a "spin liquid", namely a Mott insulator without mag- 
netic long range order at zero temperature. 

In this line of thought one would like to understand 
spin liquids, metals and their relation to superconduct- 
ing states, within one same framework. Many insights 
into this problem were obtained using mean field meth- 
ods such as slave boson mean field theory and the large-N 
expansion. [ ] More recently the development of Dynam- 
ical Mean Field Theory (DMFT) and its cluster exten- 
sions [3, 4, 5, 6] have improved the description of the 
unusual electronic structure of strongly correlated mate- 
rials. Indeed this approach applied to simple models ac- 
counts for many properties of both the (i-wave supercon- 
ducting state and the normal state of copper oxides [7, 8] 
and Ai:-organics. [ ] 

In this letter, we make a general mean-field ansatz to 
study symmetry-breaking solutions describing the insu- 
lating phases of the half-filled ti — t2 Hubbard model 
at large [/, using cluster DMFT on a 2 x 2 plaquette. 
We identify several Mott insulating phases and we take 
advantage of the substantial simplifications provided by 
the CDMFT method to obtain a simple understanding 
of their local properties. Such a program was carried out 
for the anisotropic triangular lattice in Ref. , while in 
Refs. iU, ii DMFT methods were applied to the Mott in- 
sulating regime on the kagome lattice. Unlike the above 
studies, however, our DMFT ansatz is compatible with 
dimerization, a common instability in frustrated spin sys- 
tems. 

We start with a frustrated Hubbard model on the 
square lattice 

H = - ^ (ti^-4^c^-^ + /i.c.) + /7^ni|ni|-/i^ni^ (1) 

(i,j),cr i i,cr 

where Cia- (resp. cj^) destroys (resp. creates) an electron 



at site i with spin projection a on the z-axis, tij = ti 
(resp. ^2) for pairs (i, j) of NN (resp. NNN), and tij = 
otherwise. U is the onsite Coulomb repulsion and /i is 
the chemical potential. 

We investigate the insulating phases of (1) at half- 
filling using Cellular DMFT on a 2 x 2 plaquette. To 
be more specific, we consider the following Anderson Im- 
purity Model 

^AIM = + ^ {^af3,aCiaa^(3a ^ h.C.) 

(a,/3),a 

+ {Via,acla^a + h.C.) (2) 

where 7Y|n is the restriction of the original Hamilto- 
nian (1) to one plaquette □ (the impurity), and the last 
two terms read the kinetic energy of the conduction bath 
and the hybridization between the bath and the plaque- 
tte, respectively. In substance, solving the self-consistent 
Cellular DMFT equations amounts to finding a set of 
bath parameters {e^p^a^ Via,a}a,(3,a that mimic the infiu- 
ence of the infinite square lattice surrounding plaquette 
□ in the original model (1). [ ] 

On the technical side, the use of Exact Diagonalization 
to obtain the ground state of (2) severely constrains the 
size of the bath (8 sites in the present study). Further, 
numerical tractability requires that the total number of 
bath parameters be even more limited, usually by im- 
posing additional symmetry relations among them. This 
step requires extra care, especially when investigating 
possible symmetry-breaking instabilities. Indeed, there 
is no spontaneous symmetry breaking in a finite system, 
hence the symmetry of the ground state of (2) is merely 
that of the Hamiltonian itself. Finally, attention must 
also be paid to the CDMFT self-consistency condition. 
The latter describes the embedding of the impurity pla- 
quette □ in the infinite lattice, which must be compat- 
ible with the symmetry of the thermodynamic ground 
state. These points will be discussed at greater length 
elsewhere. [ ] 

We first present the mean-field ground states of (1) at 
half-filling and for large U : in the following ti = 1 and 
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U = 16. We have checked that this U puts the system 
weh within the Mott insulating phase. In this regime the 
cost of double occupancy is avoided through the coher- 
ent hopping of two or more electrons, and to lowest order 
in ^1,2/^, (1) reduces to the Ji — J2 Heisenberg model, 
with' = 4tf//7(l + 0{{t/U)'^)). [ ] Hence the meanig- 
ful frustration parameter is (^2/^1)^. For finite [/, higher 
order processes also contribute and generate longer ex- 
change loops, starting with four-spins exchange at order 

tj/u'. 

Figure 1 shows the evolution of three order pa- 
rameters on the impurity upon increasing the frustra- 
tion ratio ^2/^1 • First, the usual Neel order is de- 
tected by computing the normalized staggered magne- 
tization (mfcj) = where ki = (7r,7r) 
and Nc = 4:. Another Neel order is anticipated at large 
frustration from previous studies of the closely-related 
Ji — J2 model (see Ref. and references therein). It 
corresponds to the entropic selection of one of the two 
degenerate impurity momenta ki = (0, tt) and (tt, 0) . Fi- 
nally we compute (D) = {S2 ■ Ss — Si ■ S2) which evi- 
dences the X — y symmetry breaking of the spin corre- 
lations. Note that {D) is non-zero if the impurity pla- 
quette carries singlets on parallel NN bonds, but also 
for the (0, tt) and (tt, 0) Neel orders. Hence the dimer- 
ization itself is detected through both i) (D) 7^ and ii) 

(^(0,7r)) (^(7r,0)) = 0. 




FIG. 1: (Color online) Evolution of the Neel and dimer or- 
der parameters computed on the impurity plaquette for in- 
creasing frustration 12/ ti (see text). Cartoon pictures of the 
various phases on the impurity plaquette are shown: in both 
Neel phases the relative orientation of the spins is indicated, 
while in the paramagnetic phases, triplets and singlets are 
pictured with parallel lines and green ellipses, respectively. 

At ^2 = we recover the (tt, tt) Neel order with a stag- 
gered magnetization renormalized to about 86% of its 
classical value. Upon increasing the frustration, {m(^^^^^) 
is increasingly renormalized until quantum fluctuations 
completely wipe out the Neel order at (^2/^1)^ 0.46. 

For 0.46 < (t2/ti)2 < 0.87 the ground state is a sin- 



glet formed by the four-spins of the impurity. Since it 
doesn't break the point symmetry of a square, we will 
refer to it as a fully symmetric paramagnet in the follow- 
ing. For 0.87 < {t2/tiY < 0.93 the ground state is stih 
paramagnetic (all Neel order parameters are zero) but it 
spontaneously breaks the rotational symmetry of the lat- 
tice and dimerizes in an almost pure product of singlets 
on the plaquette, as evidenced by I) 3/4. 

Finally SU(2) spontaneously breaks again for 
(^2Ai)^^0.93 where the (0, tt) Neel phase kicks in 
through a flrst-order transition, as suggested by the 
abrupt increase of (m(o,7r)) jumping from to about 90% 
of its classical value. This is further supported by the 
simultaneous existence of two solutions to the CDMFT 
equations in this region, namely the fully symmetric 
paramagnet and the (0, tt) Neel order. [12] 

Interpretation Cellular DMFT provides a local picture 
of each of the phases of model (1), in terms of the states 
of the self-consisent impurity model (2). 

In semi-classical phases, such as the (tt, tt) Neel phase 
at small ^2, tracing out the bath degrees of freedom re- 
sults in a reduced density matrix for the impurity plaque- 
tte that has most of its weight in one state of the plaque- 
tte. This state resembles a classical, antiferromagnetic, 
spin conflguration dressed by quantum fluctuations. The 
orientation of the spin is picked up by the self-consistent 
Weiss fleld. 

To gain more insight into the fully symmetric 
phase that takes over at intermediate frustration 
0-46 < (^2/^1)^ < 0.87, we investigate the nature of the 
two lowest eigenstates of the impurity model, noted |0) 
and From each one of these states we construct 

the reduced "density matrix" of the impurity, deflned by 
pf = TrB|Q^)(Q^| where a = 0, 1 and Trp denotes the par- 
tial trace over the bath degrees of freedom. The spectral 
decomposition of the operators pf reveals that most of 
the weight (~ 95%) is carried by a single state of the pla- 
quette noted ie. pi — \a)i ■ i{a\. This means that 
|a)i|a)B, i.e. the impurity and the bath are only 
weakly entangled in the subspace spanned by |0) and 
at stark contrast with the situation in single-site DMFT. 

Further, the nature of the impurity states \a)i provide 
a simple mean-fleld picture of the paramagnetic phases 
and their evolution with ^2/^1. Namely, in the fully sym- 
metric phase, 1 0)1 is the singlet obtained by adding two 
triplets along the diagonals (1,3) and (2,4) of the pla- 
quette, noted 1^), only lightly dressed by quantum fluc- 
tuations. Similarly, |l)i is the product of two diagonal 
singlets, noted |n). 

We will show that the disappearance of the fully sym- 
metric paramagnet at (^2/^1)^ — 0.87 is driven by the 
closing of the gap between these two singlets. We be- 
lieve that this is the local mean-fleld description of the 
abundance of low-lying singlets commonly observed at 
Quantum Phases Transition points in frustrated quan- 
tum magnets. [ ] This is supported by the observation 
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that the two impurity singlets are orthogonal (H|^) = 0, 
and span exactly the two-dimensional subspace of the 
half-filled impurity with total spin 6'i = 0. Hence, within 
plaquette DMFT there is simply no additional, linearly 
independant, half-filled impurity singlet susceptible to 
condense. 

When the singlet gap closes, the zero-energy sub- 
space is enlarged to Increasing the frustra- 
tion in 0.87 < (^2/^1)^ ^0.93 hfts this degeneracy and 
selects a particular direction within this subspace, either 
|C) = V3/2\m) + l/2p), or || I) = V^/2\m) - 1/2|K). 

Dynamical singlet correlations More insight into the 
paramagnetic phases can be gained by considering 
higher-order spin correlations. Indeed, in a seminal 
paper [ ] Wen discussed a particular class of spin 
liquids, breaking both Parity (P) and Time rever- 
sal (T), whose order parameter is the scalar chiral- 
ity {Eijk) = ('S'^ • Sj X Sk)' The relevance of these chi- 
ral spin liquids for the Ji — J2 model was claimed by 
the author, on the premise that the coherence in the 
collective dynamics of three spins on a plaquette is 
generally enhanced by a non-zero J2, possibly to the 
point where {Eij^) 7^ 0. Beyond the original mean- 
field study by Wen, an early numerical work of the 
Ji — J2 model evidenced the softening of chiral fiuctua- 
tions, although no long-range chiral order was found. [1 5] 
Within the CDMFT treatment of the related Hubbard 
model (1), the weak entanglement between the bath 
and the impurity yields {{)\Eijk\G) 0^ {M\Eijk\M) in the 
fully symmetric phase. Consider for instance the triplet 
(1,2,3) of impurity sites. Direct computation shows that 
^123 |K) = iVS/Ma) and £;i23p) = -iV^/MM). Hence 
(0|£^i23|0) (X = in the fully symmetric phase 

and a similar computation in the dimer phase yields 
(0|£^i23|0) (Z|£^i23|Z) = as weh. Thus there is no 
instability towards chiral symmetry breaking in our ap- 
proach either, and we need to address the dynamics of 
£'123. Namely we compute the following dynamical cor- 
relation 

(£123(^^)^123) = Yl M^i23\0)\^5{u; - e, + eo) (3) 
q 

where the sum runs over all eigenpairs {^q^\q)}q of the 
impurity model (2). The above discussion shows that in 
the fully symmetric phase, the lowest pole in Eqn. 3 is 
located at the singlet gap energy uoi = ei — eo with weight 

Zi = Kl|i?i23|0)|2c.(\/3/4)2|B(l|0)Bp. 

In Figure 2 we show the evolution of coi and Zi with 
increasing ^2/^1 . As expected, increasing the frustration 
^2/^1 from ^2 = reduces uoi and increases Zi, evidenc- 
ing the enhancement of coherence of the three spins on 
a plaquette. In both paramagnetic phases the weight Zi 
decreases very slowly, in the 0.16 — 0.155 range, showing 
that |b(1|0)b| — 0-9- Hence, not only are the bath and 
the impurity weakly entangled, but also the bath states 
that accompany the lowest lying impurity states are quite 
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FIG. 2: (color online) Evolution of the frequency ui and in- 
tegrated weight Zi of the lowest pole of the chiral correlation 
{Ei23{oo)Ei23) with increasing frustration. 



similar. In the following, this will allow us to ignore the 
bath components of |0) and |1) temporarily. Although 
there is no significant "decoherence" induced by the cou- 
pling of the impurity to the bath, the bath does play a 
crucial role in determining the phase diagram: without 
it we would not be describing a system in the thermo- 
dynamic limit, able to sustain broken symmetry phases 
and transitions between them. 

The upper panel of Figure 2 shows that for 
0-46 < {t2/tif < 0.87, the singlet gap ui between |^) 
and !□) decrease linearly towards zero as (^2/^1)^ 0.87. 
For comparison, the lowest pole of the cluster dynami- 
cal spin-spin correlations {mk^iuo) • rri-kj) is obtained for 
ki = (tt, tt) in this region and yields a local estimate of 
the spin gap of about 0.2 ti. 

In the fully symmetric phase, any combination 
l"^) = a\M) -\- b\\Z\) is not an eigenstate and instead 
has a dynamics with characteristic frequency coi. As 
(^2/^1)^ 0.87, a;i ^ and all such states become eigen- 
states with zero-energy. This obviously includes chi- 
ral states, defined by (^|£^i23|^) 7^ 0, and such that 
6a* — a6* 7^ 0, but also the above-mentioned dimer states 
1 1 I) and one of which will be eventually selected upon 
increasing ^2/^1 above 0.87. 

We now ask whether these chiral and dimer states 
exhaust the possibilities of phases that can be de- 
scribed by plaquette DMFT and which compete near 
(^2/^1)^ — 0.87. One route consists in extending the 
above discussion of the chiral correlations. Namely we 
are looking for any additional hermitian operator A 
such that the correlation {A{uj)A) has its lowest pole 
at the energy of the singlet gap coi in the fully sym- 
metric phase. According to Eqn. 3 this translates as 
{0\A\0) {m\A\m) = and {1\A\0) {n\A\m) + 0. In 
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such a reduced subspace the possibihties are very hmited: 
the set of 2 X 2, hermitian, off-diagonal matrices is exactly 
spanned by the two Pauli matrices (Jx^cFy In particular, 
the above computation of £^123 can be recast in opera- 
tor form as £'123 = a/3/4 cfy. A similar computation for 
the dimerization operator D leads to D = a/3/2 cFx- In 
this synthetic form it is clear that we have exhausted the 
list of linearly independant operators whose fluctuations 
become infinitely soft at the transition point. 

Interestingly, the local mean-field picture developed 
above allows us to further elaborate on the nature 
of the chiral phase that competes with the dimerized 
phase near (^2/^1)^ — 0.87. Indeed, repeating the above 
computation for all triplets of impurity sites leads to 
^123 V3/4 cFy = -£^234 = £^341 = -£^412- It is inter- 
esting to note the chiral spin liquid proposed by Wen 
in its large- treatment of the Ji — J2 model obeys the 
same spatial symmetries. [ ] Further, Wen derived both 
i) a chiral liquid and ii) a dimer solution of the mean- 
field equations. Both solutions were found to compete 
at large frustration J2 — Ji^ although dimers are always 
energetically favored, in qualitative agreement with the 
present study. 

Conclusion It is remarkable that a mean-field theory 
based on a plaquette is able to describe the phases of 
the Ji — J2 model, as obtained within more exact ap- 
proaches. [13, 16] The present study provides simple car- 
icatures of the local aspects of these phases. Mean- 
field theory, however, does not contain the effects of 
long-wavelength fluctuations and is not expected to de- 
scribe the location of the phase boundaries accurately. 
CDMFT also ignores more complicated forms of long- 
range order which cannot be tiled with impurities. Over- 
all we expect that, upon increasing the size of the 
impurity, the phase boundaries in Figure 1 will shift, 
moving the paramagnetic region 0.46 < (^2/^1)^ < 0.93 
closer to the values of the nearby Ji — J2 model, namely 
0.4 < J2/J1 < 0.6. [ "] [1^^]. Nevertheless, we expect the 
local physics described in this paper to be robust and play 
an important role both at finite temperature and finite 
doping, which deserves further study within plaquette 
DMFT. In particular, doping this system in the region 
with low-lying chiral fluctuations may generate orbital 
currents or a d-density wave. Although no instability to- 
wards either of those orders was detected in the ^2 = 
case at finite doping, [ ] the present study shows that 
this issue should be reconsidered for ^2 7^ 0. 

On the experimental side, the search for physical re- 
alizations of a ti — t2 one-band Hubbard model resulted 
in materials with either small ^2/^1, such as cuprates, 
which have (7r,7r) Neel-order at low-temperature, or 
very large ^2/^1, such as Li2VOSi04, lying well in the 
(0, tt) Neel phase. [ , xu] However, the recently syn- 
thetized material PbV03 [20] has intermediate frustra- 
tion [21] and offers a possible realization of the non- 
magnetic phase described in this paper. In particular. 



Raman scattering under pressure in various geometries 
should be able to probe for the low-lying chiral fluctua- 
tions [ , ] and their evolution with frustration. 
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